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BEC AND THE NEW 
WORLD OF COHERENT 
MATTER WAVES 

Allan Griffin 

Department of Physics, University of Toronto, 
Toronto, Ontario, Canada M5S 1A7 

After decades of effort to produce an atomic Bose condensate, this was 
finaUy achieved in June, 1995 in a laser-cooled magnetically trapped gas 
of *^Rb atoms by Eric Cornell and Carl Wieman at JILA (University of 
Colorado and N.I.S.T.). As of July, 1999 (4 years later), there are now 
over 20 experimental groups around the world who can routinely produce 
and study such atomic Bose condensates. In addition, over 1000 theoretical 
papers have been published! 

I will give three introductory lectures on this new Quantum Phase of 
Matter, which I think will continue to be a growth point of fundamental 
physics research for the next decade and will also be the source of new 
technologies based on using this source of coherent matter waves. A brief 
sketch of my three lectures is as follows: 

1. A BRIEF HISTORY OF BEC STUDIES. 

• Before 1995, going back to pioneering work of Einstein (1925) 
and Fritz London (1938). 

• Introduce the concept of a Bose macroscopic wavefunction $(r, t) 
describing a Bose condensate and the key Gross-Pitaevskii equa- 
tion of motion for ^(r, t). 

• Review some of the results obtained over the last 4 years - and 
indicate why trapped Bose gases are so interesting, possibly even 
more than superfluid "^He and BCS superconductors, which also 
involve Bose condensation. 

2. DYNAMICS OF A PURE CONDENSATE (T < Tbec)- 

• Time-dependent Gross-Pitaevskii (1961) equation of motion. 

• Crucial effect of weak interatomic interactions. 

• Collective oscillations of the condensate. 



• The Bose condensate as a classical quantum object! 
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3. DYNAMICS OF COUPLED CONDENSATE AND NON-CONDEN- 
SATE COMPONENTS (TWO-FLUID HYDRODYNAMICS): 

• Derivation of a quantum Boltzniann equation for the non-conden- 
sate excited atoms and a generalized GP equation for the con- 
densate. 

• More complex behaviour than in superfluid *He. 

• Comparison with the well-known Landau two-fluid theory (1941). 

Further references which are relevant for these lectures are: 

1. BEC Homepage (maintained by a BE C theorist, Mark Edwards) 



bttp : / /amo . phy. gasou . cdu/bec . ht ml/ . This widely- used homepage has 



played a crucial role in BEC research. 

A recent article [|l| by Dalfovo, Giorgini, Pitaevskii and Stringari is 
an authoritative review of recent theory on atomic Bose condensates 
and expands on the material I cover in Sections |l| and ^j. 

Articles in a book §, Bose-Einstein condensation in atomic 
gases, Proceedings of the International School of Physics "Enrico 
Fermi", ed by M. Inguscio, S. Stringari and C. Wieman. This contains 
many review articles on current BEC research. In particular, I call 
attention to (all these can be downloaded from the LANL website 
under cond-mat): 

• W. Ketterle, D.S. Durfee and D.M. Stamper-Kurn, "Making, 
probing and understanding Bose-Einstein condensates" - a 100 
page review of recent experiments . 

• A. Griflin "A brief history of our understanding of BEC: From 
Bose to Beliaev" 0. 

• A. Griffin, "Theory of excitations of the condensate and non- 
condensate at finite temperature" [Q. 

• A.L. Fetter, "Theory of a dilute low-temperature trapped Bose 
condensate." A very detailed analysis at T = 0. 

A long article by Zaremba, Nikuni and Griffin |^ on the non-equilibrium 
dynamics of trapped Bose gases at finite temperatures. This expands 
on the material covered in Section ^ of these lectures. 



1 AN OVERVIEW OF PAST AND RECENT 
WORK 
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1.1 Some history before 1980 

Einstein predicted that a non-interacting gas of atoms (Bosons) would un- 
dergo a phase transition at low temperatures, when a macroscopic {Q{N)) 
number of atoms occupy the lowest energy level (in a uniform ideal Bose 
gas, this is the the zero momentum single-particle state ) . His work was in- 
spired by a novel derivation of the Planck distribution for photons by Bose 
in 1924. The basic physics of this phase transition is worked out in every 
text in statistical mechanics The simplest way of estimating Tbec is 
to note that the transition occurs when 

d — average distance between atoms ~ . , (1-1) 

where At is the thermal De Broglie wavelength of a gas of atoms at tem- 
perature T, 

\mkBT J 



The criterion in (1.1) is equivalent to nX^ > 1, while a more careful analysis 
[I gives nX^ = C(3/2) = 2.612. One sees that At ^ large as T ^ 0. When 
At ^ d, all the atoms become correlated and the gas exh ibits new collective 



behaviour (even in absence of interactions). Using (IT), one finds that 

kBTBEC-'^-^n'"- (1.3) 
m 

Below Tbec, the number of atoms in the p = single-particle state in- 
creases and is given by the well-known formula 



Nc{T) 
N 



J, \ 3/2' 



Tb 



EC 



(1.4) 



At T = 0, all the atoms in an ideal gas are in this p = state (this state 
is the Bose condensate in non- interacting 3D gas). 

Nothing much happened until 1938. Then the neglected work of Einstein 
was re-discovered and developed by Fritz London |^], who suggested that 
it might be the basis of an explanation for the strange effects noticed in 
liquid ''He at ~ 2.17i^. London's suggestion was based on the fact 
that the ''He atom was a "composite" Boson {S = 0) and the formula in 



( |l.3| ) gives Tbec ~ 3iir if we used the density for liquid ''He. L. Tisza 
used London's idea and suggested (somehow!) that the condensate atoms 
act in a coherent way - a new collective degree of freedom moving without 
friction. This "picture" led to a rudimentary two-fluid model that could 
explain experiments showing superfluidity (especially by Kapitza as well as 
by Allen and Meisner) as a counterflow of the superfluid and the normal 
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fluid. In a dilute, weakly interacting Bose gas, these two components are 
the condensate and non-condensate, respectively. 

Until the 1960's, the theory of BEC in interacting systems was dominated 
by efforts to use it to understand superfluid ^He. The London-Tisza scenario 
was essentially correct, but to formulate it properly needed field-theoretic 
many body techniques and the concept of broken-symmetry, which were 
only developed in the period 1957-1965. In this period, a large amount of 
work was done on a toy problem, a dilute weakly interacting Bose gas, since 
a Bose liquid like superfluid ''He was too difficult to deal with theoretically. 
These early studies are the foundation of our current understanding of 
trapped atomic gases. 

A phenomenological theory of superfluid *He was introduced by Landau 
in 1941, based on the idea of quasiparticles (phonon-roton spectrum) and 
a two- fluid superfluid hydrodynamics This brilliant theory has been 
very successful and is the basis of modern descriptions of superfluid "^He 
[Eo|. However, it made no explicit mention of BEC or even the fact that 
^le atoms obeyed Bose statistics. Only in the 1960's did it become clear 
that Landau formulation had it's microscopic basis in the existence of a 
condensate macroscopic wavefunction, 

<i>(r, t) =< ^(r) v'iV:(M)e''^'"'*\ (1-5) 

where ip{r) is the quantum held operator (see Section ||). This concept was 
first formally introduced by Beliaev in 1957 ||ll|, extending the pioneering 
work of Bogoliubov in 1947 ||l2|. The superfluid motion is associated with 
the gradient of the phase of this two-component order parameter, 

^ ^^{0o+r■Ve) . = V0 = ^ (1.6) 

Ti 

The essential relation between superfluidity and $(r,t) is simply and ele- 
gantly described in the classic monograph by Nozieres and Pines ||l^ . For 
further discussion of the development of our current understanding of Bose 
condensates, see the review article by Griffin in Rcf. ||^. 

It might be useful to make a brief digression here on the BCS theory of 
superconductors based on formation of Cooper pairs (total spin 5 = 0), 
which was developed in 1957. In the BCS theory, superconductors exhibit 
the same kind of macroscopic quantum behaviour as superfiuid *He, as first 
argued in the late 1930's by F.London. In the Gor'kov version of the BCS 
theory, the spin-singiet Cooper pair order parameter 

<i>(r,i) =< ViT(i-)V'i(r) > (1.7) 

is the equivalent of 

$(r, i) =< ^/.(r) > (1.8) 

in Bose superfluids. The essential equivalence of these two systems was ob- 
scured by the complexity of the original many-particle BCS-wavefunction, 
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due to the large spatial size of overlapping Cooper pairs. As a result, Cooper 
pairs form and condense into a coherent state at same temperature. It was 
only in the 1980's that theorists (Leggett, Nozieres and others jlj]) realized 
that for small Cooper pairs (tightly bound and hence high Tc), the BCS 
theory smoothly goes over to theory of a weakly interacting Bose gas of 
non-overlapping Cooper pairs. Recently a gas of '^"K atoms (a composite 
fermion) have been laser-cooled at JILA to slightly below the Fermi tem- 
perature jl^ . There are ways to make the interaction attractive (see later) 
and then one might look for the formation of Cooper pairs |l6| of ultra-cold 
Fermi atoms! 

1.2 More recent developments (1980-1995) 

Since the 1970's, there has been increasing interest by experimentalists to 
find a "pure" form of BEC, namely, in a low temperature gas. The two early 
candidates for Bosons were excitons (electron-hole pairs) in semiconductors 
and spin-polarized hydrogen atoms (see articles by Greytak and by Wolfe 
et al. in the book mentioned in Ref. [0). A gas of atoms was predicted 
to be stable as a gas even at T = 0. This is because the atoms cannot 
combine since there is no bound state of the interatomic potential between 
two spin-polarized H atoms. Thus, one cannot form liquid or solid phase. 
Many people got interested in BEC in iJ| gas, including theorists p7[ |. 

Several of the key ideas that led to success in alkali atoms in 1995 grew 
out of the pioneering work on 7J| gas in the 1980's. However 3-body inter- 
actions become increasingly important at higher densities and these cause 
spin flips, allowing formation of H- molecules. High densities were needed 
since cooling was by cryogenic methods, which could reach ~ 10~'^K but 
no lower. BEC in iJ| gas was finally produced at MIT in June, 1998, after 
almost 20 years of work jl^. Unfortunately, the alkali atom condensates 
are much easier to create and study, and appear to be more interesting 
gases. 

Since the early 1990's, attention has focussed on the alkali atoms: Li, Na, 
K, Rb, Cs. The atoms are Bosons, with an even number of neutrons. The 
strategy was to use laser-cooling to get to very low temperatures, where 
the low density gas would Bose-condense. The essential idea behind laser 
cooling is that when an atom absorbs a photon, it slows down. In the 
summer of 1995, BEC was announced by three groups led by: 

• C. Wieman and E. CorneU (JILA), using ^''Rb atoms [|l9|. 

• W. Ketterle (MIT), using ^SNa atoms 

• R. Hulet (Rice University), using ^Li atoms ||2l| . 

Parenthetically, it is now the general feeling that the original Rice data, 
interesting as it was, did not give an unambiguous signature of a (very 
small, since the interaction is attractive) condensate in ^Li gas p^ . 
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Alkali atoms are perfect for BEC studies. They have a magnetic moment, 

and hence can be trapped by magnetic fields. They essentially have a "one- 
electron" structure. They are thus simple atoms, and have been well studied 
by atomic physicists. One can easily selectively flip the "spin" of higher 
energy trapped atoms. These "hot" atoms are then quickly ejected from 
the magnetic trap and the remaining atoms quickly thermalize to a lower 
temperature. This "evaporative cooling" is very efficient and quickly brings 
one into the temperature region required for BEC. 

It is useful to mention a few experimental facts about the magnetic traps 
currently in use. As it turns out, these traps are well described as a har- 
monic potential 

Vex{r) = ]^muolr^ (isotropic) (1.9) 

= \m{i^lxX^ + ^lyV^ + '^Izz'^) (anisotropic). (1.10) 

Most current traps are either: 

pancakes , uJqz 3> i^ox{~ ^oy) 

cigars , lOoz < i^ox{= ^oy) (1-11) 

and the trap frequencies are of the order t^o ^ 27r x lOOHz. In 1995, the 
first condensates were small ~ 10^ atoms and Tbec ~ lOOnK. However 
in 1999, the condensates can be quite large ^ 10^ atoms at Tbec ^ M^. 
These have a size ~ many microns, which can be easily seen optically . 
When the condensates are small, the trap is turned off and cloud allowed 
to expand, and then measured by optical methods. The results are simple 
to analyze if gas is non-interacting. However, more analysis is needed to 
include the effects of interactions during expansion. 

Early reports on atomic condensates discussed the system as an ideal 
Bose gas. It was soon realized that even in these very dilute gases, the 
interactions played a crucial role. Indeed, since the 1960's, it has been 
understood that an interacting Bose gas is quite different from an ideal 
Bose gas. In particular, interactions stabilize (or "lock") the phase of the 
condensate and allow coherent properties to emerge. (We recall that a free 
Bose gas has a condensate but is not a superfluid). Of course, interactions 
are also crucial for cooling. After hot atoms arc removed by rf-induced 
spin-flips, it is important that remaining atoms can quickly re-thermalize 
through collisions. 

However, it is useful to first consider an ideal Bose gas in a trap, to 
illustrate some characteristic features. For atoms in an external potential, 
we have 

N = Y,h{^i\ (1-12) 
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where the Bose distribution is 



In a harmonic trap, the energy levels are: 

3 

ei = en,ri„n, = {nx +ny + n^ + -)huJo, (1-14) 

with Tlx, ny^jiz = 0, 1, 2, 3 . . .. The condensate is described by the ground 
state single-particle wavefunction 

</>o(r)~e-'-'/2--o, (1.15) 

where the S.H. oscillator length is ano = i^T^jj^ ^ IfJ^m in current 

traps. Clearly ano gives the "size" of condensate nco(r) — |<i>o(r)P in a 
trapped gas. For the non-condensate density no(r) (often called the "ther- 
mal cloud"), we can use the semiclassical limit, since the thermal energy 
[ksT) is much larger than the spacing between the S.H. energy levels {fiuJo). 
Then we have 



e-'-'/^fl^ (1.16) 



where 



= \ 2 = "-ffo t — ~> ano- (l-l') 



^0 



Thus we see that the size of the thermal cloud Rt is much larger than 
the condensate. The signature for condensate is this sharp high density 
peak at the centre of the trap, which suddenly starts to grow out of the 
broad thermal distribution at the predicted transition temperature Tbec- 
As T ^ (effectively T < 0.4 Tbec), the thermal cloud stead ily disap- 



pears as all atoms go into the ground state (?!)o(r) given by ( [1.15[ ), which is 
the macroscopic wavefunction for a non-interacting trapped Bose gas. The 
temperature of the gas is measured from the temperature dependence of 



the tail of the thermal distribution given by (1.16) 



It is easy to calculate the transition temperature for atoms in an har- 



monic trap. Separating out the condensate contribution in (1.15), we have 



iV = iVc + ^/°(e»); iVc= / dr|0o(r)p. (1.18) 

Making a change of variable ptiujonx = fix, and using the continuum ap- 
proximation, we have 

/k /"OO /"OO /"OO 

N — Nr ~ I / dfiT / dfioi / dn 



fi^^a J Jo Jo Jo e("-+"!'+"-) - 1 
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C(3) 



(1.19) 



We note that the chemical potential /ig = jfiujQ^ but ^ 1 and hence 
the zero point energy has been neglected in ( 1.19| ). Since A^c = at Tbec-, 
we have 



B-lBEC 



0M{N3)hujQ, with 



Nc{T) 
N 



\Tb 



EC 



(1.20) 



One can improve on these simple estimates for Tbec and Nc{T), but the 
corrections are only a few percent at best . 

Interactions make a dilute, weakly-interacting Bose condensed gas into a 
full non-trivial many body problem, even though the system is very dilute 
and the interactions are weak. In a dilute gas, we need only consider binary 
collisions. The real interatomic potential v{r) has a hard core with a radius 
of a few Angstroms and a weak, long-range attractive tail. In a dilute, 
very cold gas, we can approximate v{r) using the s-wave scattering length 
approximation effectively replacing v(r) by a pseudopotential l23] 



v{r) 



a6{r) = g5{r). 



(1.21) 



We require a ^ average distance between atoms, or na'^ <C 1, which is very 
well satisfied in these gases. For alkali atoms, v{r) almost has a bound state 
of two atoms. This quasi-bound state is very sensitive to the long range part 
of the potential, and thus the value of the s— wave scattering length a can 
be very large. Current values for atoms used in BEC experiments are: 



87 



23 



Rb 

Na 
^Li 



a = 58 A 
a = 28 A 
a = -14 A 



(1.22) 



One can adjust the energy of the quasi-bound state and, as a result, 
change the value of s-wave scattering length a with a small magnetic field. 
Near a so-called Feshbach resonance, one can even change the interaction 
sign, going from repulsive (a > 0) to attractive (a < 0). There is a lot of 
current work [ p4| trying to exploit this ability to change the interaction 
strength and sign by "turning a knob" , perhaps even making a —^ 01 

The key equation for the macroscopic wavefunction for a T = conden- 
sate was written down and discussed by Pitaevskii |^ and Gross in 
1961: 



ih 



a$(r,t) 
dt 



2m 



(1.23) 
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where Uc = I'&P- This equation describes the condensate atoms moving in 
dynamic seh-consistent Hartree field produced by the condensate, 

Vnir, i)^ j dr'v{r - v')n,{r\ t) = gn,{v, t). (1.24) 



The non-hnear GP equation ( 1.23| ) wiU be the subject of Section ||. Hun- 



dreds of papers have been written on it in the last four years. For T < 0.4 TbeC: 
it describes both the static properties and the dynamic fluctuations (linear 
and non-linear) very well, usually within a few percent jl], 

To complete this brief introduction, I mention several other research 
topics that make trapped Bose-condensed gases so exciting: 



1. Alkali atoms have several different atomic hyperflne states. Apart 
from Cs, the alkali atoms have a nuclear spin 1 = 3/2 and an electron 
spin S = 1/2. Thus the total spin operator F = I + S has values of 1 
and 2, leading to 8 possible atomic states. One usually works with one 
of these atomic hyperfine states - trapped in a magnetic well. How- 
ever, purely optical traps (MIT) using the dipole force of a laser beam 
can be used to trap low energy atoms in different hyperfine states. 
Thus one can now also deal with gases with several different states, 
ie, a multicomponent Bosc gas. Moreover, one can induce transitions 
between different hyperfine states. One sees that spin is a new degree 
of freedom in such multicomponent Bose-condensed gases ^ . 

2. Independent of the special "coherent" features of a Bose condensate, 
these trapped gases give us a source of high density, very cold atoms. 
Lene Hau |2^ has used this high density to slow down the speed of 
light to that of slow car (~ 40 km/hour) using self- induced trans- 
parency. One can also switch on an optical lattice (produced by in- 
tersecting laser beams) on a trapped Bose gas Turning off the 
magnetic trap, the low energy atoms will occupy the potential min- 
ima of this periodic lattice. This could not be done with high energy 
atoms since the dipole-induced potentials of the optical lattice are 
very weak. With ultra-cold trapped fermions, one may also be able 
to produce a Hubbard model, of the kind extensively studied in con- 
nection with the cuprate-oxide high temperature superconductors. 

We agree with the opinion of Pitaevskii that the discovery of BEC 
in alkali gases "can be considered as one of the most beautiful results of 
experimental physics in our century" . The next two sections will fiesh out 
this qualitative overview with some theoretical calculations on the collective 
oscillations at T ^ Tbec (pure condensate) and at finite temperatures 
T Tbec (mixture of condensate and non-condensate) of these strange 
quantum "wisps of matter" . 
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2 DYNAMICS OF THE PURE CONDENSATE 



The theory of interacting Bose-condensed fluids is most usefuUy discussed 
using quantum field operators. This procedure was formahzed by Behaev 
(1957) and developed by Bogoliubov |Q, Gavoret and Nozieres Martin 
and Hohenberg and others in the 1960's [Q. We recall: 

if]'^ (r) = creates atom at r 

■0(r) = destroys atom at r. (1-25) 
These fields satisfy the usual Bose commutation relations, such as 

V^(r),?/'+(r')j =(5(r-r'). (1.26) 

All observables can be written in terms of these quantum field operators, 
such as the interaction energy 

Ve^t = iy"dry"drV+(r')^A+(r)v(r-r')i/'(r')^(r) 

^ ^9 J rfri/i+(r)v3+(r)V;(r)V3(r). (1.27) 

The crucial idea due to Bogoliubov (1947) and later generalized by Be- 
liaev is to separate out the condensate part 

^(r) = (^(r))+^(r), (1.28) 

where 

■0(r)^ = $(r) = Bose macroscopic wavefunction. (1-29) 

This quantity plays the role of the order parameter for the superfluid phase 
transition: 

$(r) =0 T>Tc 

7^0 T<Tc. (1.30) 

We note that <&(r) = \fWc&^^ is a 2-component order parameter. Clearly, 
$(r) is not simply related to the many-particle wavefunction ^(ri, r2, . . . tat). 
The thermal average in < ■0(r) > involves a small symmetry-breaking per- 
turbation to allow $ to be finite, 

BsB = Jdr [7y(r)v3+(r) + 77*(r)^/>(r)] . (1.31) 

It is useful to make a few comments on the physics behind $(r, t). <i>(r, t) 
is a coherent state, with a "clamped" value of phase - rather than a Fock- 
state of fixed N, with no well-defined phase. 3>(r, t) acts like a classical 
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field, since quantum fluctuations are negligible when Nc is large. Probably 
P.W. Anderson deserves the greatest credit for understanding (in the period 
1958-1963) the new physics behind working with a broken-symmetry state 
$(r, i), both in BCS superconductors and in superfluid ''He Q. It cap- 
tures the physics of the new phase of matter (such as the occurence of the 
Josephson effect) and the associated superfluidity. The symmetry-breaking 
perturbation allows < ■0 > to be finite. More precisely, it allows the system 
to internally set up off-diagonal symmetry-breaking fields, which persist 
even when the external symmetry-breaking perturbation in (1.31) is set to 
zero at the end (77 ^ 0). The same sort of physics is behind the BCS theory 
of superconductors. 

The exact Heisenberg equation of motion for the field operator is 



ih 



dt 



2y72 



-^ + VUv)+5V{v,t) 
Zm 

+ f7(r) + 3-0^ (r, i)V'(i", i)V'(i", t), 



0(r,t) 



(1.32) 



where 6V{r,t) is a small time-dependent driving potential. This gives an 
exact equation of motion for $(r, t) = ( V'(r, t) 



ih 



a$(r, t) 
dt 



2yj2 



+ VeAr) + SV{r, t) 



2m 

+ ?7(r)+g(v3+(r,i)0(r,t)V(r,t)), 



$(r, t) 



(1.33) 



with 



= |$|2$ + 2|$p0 -I- $2^+ -1- $*Vi0 -t- 2<i>4>+'4; + (1.34) 
Taking the symmetry-breaking average, one finds 

l^jj'^ip'ip'^ = nc$ -1- m$* + 2n$ + (^'ip+ ipip^ , (1.35) 

where 

nc(r, i) E 



m(r, t) 



^{r,t)f = condensate density 
?/'^(r, t)-0(r, <)^ = non-condensate density 

^(r, i)'i/j(r, <)^ = off-diagonal (anomalous) density. 

Here we have separated out the condensate and non- condensate parts 

7/, ^ + = * + V'- (1-36) 



In general, the equation (1.33) for $(r, i) is not closed - it is coupled to 
the dynamics of the non-condensate. However, in this Section we limit 
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ourselves to T ^ Tbec, where we can assume the non-condensate fraction 
is negligible, leaving 



ih 



9$(r,t) 
dt 



2y72 



2m 



+ K.(r)+g|<i>(r,i)|2 



*(r,i). 



(1.37) 



This is the famous time-dependent Gross-Pitaevskii equation for the con- 
densate macroscopic wavefunction. It gives a complete description of the 
dynamics of a coherent matter wave at T = 0. 



2.1 Static condensate 



We first consider the time-dependent stationary GP equation, which has 
the solution 

(^(r, t)) = $(r, t) = $o(r)e-^^*/^ (1.38) 
where ^ is the chemical potential. The physics behind this can be seen from 

(^N - l\^{r,t)\N'^ ^ e*^"-i*/''(|7V-l|i/'(r)|iV^e-'-^"*/^ 



Ne' 



Using (L3|) in (1.37) gives 



2y72 



rv 

2m 



+ V,,ir)+g\Mr)\' 



*o(r). 



(1.39) 



(1.40) 



The static GP equation for the static condensate wavefunction ^o{r) is 
thus 



$o(r) = 0. 



(1.41) 



A simple approximation in solving (1.41) is to ignore the kinetic energy 
of the condensate, ie, neglect the — term. This is called the "Thomas- 
Fermi" approximation (TF) in the recent Bose gas literature M . In this TF 
approximation, the static GP (1.41) equation for ^o{r) reduces to 



[ye.(r)+5l'i>o(r)P] 
which is easly inverted to give the condensate density profile 

1 



(1.42) 



- [/U - Vexir 
9 

1 



2 2 



> 0. 



(1.43) 
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Clearly in the TF approximation, the "size" of the condensate is Rtf, 
where 

M = T^rnujlR^p. (1.44) 



One finds /i from the condition J (irnc(r) = Nc = N, or 



An 



drr^ — 
9 



This gives 



15- 



Na 



2/5 



a-HO 



(?i/mwo)i/2. 



(1.45) 



(1.46) 



We note that SH osciUator length gho is the size of the ground state 
wavefunction ( 1.15| ) of an atom in a parabolic potential. Combining (1.45) 
and (1.44) gives 



f Na\ 

Rtf = aHo 15 

V aHoJ 

Na 

> anoi II > 1- 



(1.47) 



We thus find the surprising result that interactions (while weak) spread 
out the ideal gas condensate {Rtf ^ a^o) and decrease the density of the 
condensate at centre of trap. The TF approximation for nco(r) is very good 
for large TV, except for a small region near the edge of condensate (~ Rtf)- 
Experimental data confirms these CP predictions for the nco(r) condensate 
profile, emphasizing that the condensate is not simply the ground state 
wavefunction of the harmonic trap potential. 



2.2 Dynamics of the condensate (collective modes) 

If we linearize around the static equilibrium value of the condensate 

$(r, t) = e-'^*/'* [$o(r) + S^r, t)] , (1.48) 

where ^$ <C <i>o, we see that 



ih 



9$ 



dt 

which gives 
.5(5$(r, t) 



2m 



+ Ve^r) + g [|$o|' + + *o<5$*] [$o ■ 



e K 



(1.49) 



ih- 



dt 



2m 



+ Kxt(r) + 25|$o|'-M 



(1.50) 
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We also have a similar equation of motion for S^*{r,t). Solving these two 
coupled equations with the ansatz 



(5$(r,i) = u(r)e-'"* + w(r)e'' 



(1.51) 



we find two coupled "Bogoliubov equations" for the amplitudes u and v 



+ Vexir) - /i + 2gnco{r) u{r) + gnco{r)v{r) = E,u{r) 



+ Vex{r) - Ai + 25nco(r) v{r) + gncoir)u{r) = -E,v{rXl.52) 
Here Ei = Twj are the excitation energies of the condensate. 



The equations in (1.52) have been solved numerically by several groups 
and the observed oscillations are in good agreement with these predictions 



|3^. As an illustration of the physics, it is useful to solve (1.52) for a 
uniform Bose gas. In this case we have 



u{y) = ue'''''' 
f (r) = ve'^-^ 



(1.53) 



which gives 



{hujf 



2 1,2 



2m 



— fi + 2gn 



cO 



cO 



el + 2gncQek. 



(1.54) 



This is the famous Bogoliubov spectrum at T = Here we have 

used fiQ — grico discussed earlier (see (1.42)). One finds a phonon region at 
long wavelengths 



V TO / 



1/2 



(1.55) 



The cross-over from particle-like to this collective phonon region occurs at 
kc, where 

7j2^2 

— — ^ = 2gnco ^ kc = y/Armicog/h. (1.56) 
2to 

This shows how the interactions changes the qualitative nature of low en- 
ergy exciations in a Bose-condensed gas. This feature can be shown to 
stabilize superfluid motion against dissipation [ p^ . 

These oscillations of the condensat e can be understood as excitations 
involving the non-condensate. Using ( |l.36 ), the Hamiltonian is given by 



H~fiN 



2y72 



rv 

2m 



+ Ka,-(r) -/i 



Vi(r) 
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2y72 



2m 

2.9 j dr|$o(r)PV>(r)^(r) 
/ dr$2(r)^+(r)^+(r) 



+ Vex (r) - 



*o(r) 



+ dr$*2(r)v;(r)V;(r). (1.57) 

We can diagonalize the quadratic part of the Hamiltonian, using 

^(r) = h + <(r)a+] , (1.58) 



where [dijcij^] — 6ij (Boson quasiparticles). Thus one finds 
H — fiN = const. + hujiafcti. 



(1.59) 



This transformation shows how the non-condensate part of Hamiltonian 
can be reduced to a system of non-interacting quasiparticles with a spec- 
trum identical to the condensate fluctuations. This equivalence is easy to 
understand. The condensate fluctuations 



(1.60) 



(5$= (V(r))-$o 



can be calculated to first order in the symmetry-breaking perturbation 



H.^h = dr 



rjtjj + r]*'ip 



Then standard linear response theory H] gives (schematically) 



(1.61) 



< > 77+ < [■0, tp] > ri* 



(1.62) 



This shows that the single-particle Green's functions of the non-condensate 
fields have the same spectrum as 5$. This identity of the spectrum of den- 
sity fiuctuations and single-particle excitations is a characteristic signature 
of all Bose-condensed systems which persists at finite temperatures |Q, || . 

One interesting collective oscillation is the dipole mode corresponding to 
rigid oscillation of the centre of mass of the static condensate profile, and 
predicted to have the trap frequency loq. This mode is described by 



nc{r,t) ^ nco{r - r]{t)), r]{t) ^ Vc, 



(1.63) 
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where the time-dependent centre of mass satisfies 

^--oVO- (1.64) 

This mode at frequency ujq is special feature of a parabolic trap and is called 
the Kohn mode for in the case of interacting fermions . This "sloshing 
mode" is used in BEC experiments to measure the natural frequency wq of 
the trap and it exists at finite temperatures as well (see Section ^ . 



2. 3 Quantum hydrodynamic formulation 



One often rewrites the time-dependent GP equation using the amplitude 
and phase variables ||l| 

$(r,i)=V^e^^ (1.65) 



Inserting this into the GP equation ( 1.37 ) and separating out the real and 
imaginary parts of the equation gives: 



dnc{r,t) 
dt 



V • nc(r, i)vc(r, t) — 0, continuity equation (1.66) 



/ic(r,<) + ^mvj{r, t) 



, Josephson equation(1.67) 



Here the gradient of the phase is related to the superfluid velocity by 



mvc(r,t) = hV9{Y,t) 
and the condensate chemical potential is 

/rTr 



Hcir,t) = -- 



2in^/n, 



+ Vexir) + g72c{r, t). 



Taking the gradient of ( 1.67 ) gives 



m 



dvc 

at 



-Vv; 



(1.68) 



(1.69) 



(1.70) 



The equations in (1.6G) and ( 1.70| ) "look" like those in classical hydrody- 
namic theories. They show that the condensate can be described in terms 
of coherent motions involving two variables: 



nc(r,t), Vc(r,t) 



(1.71) 



The Landau 2-fluid equations |9[ ^ reduce to these same equations at 
T = (where ps = p, Pn = 0), namely 



— + V • nvc = 

ot 



1. 



(1.72) 
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We will find these equations useful in Section ^, where they complement 
the hydrodynamic equations describing the non- condensate in the collision- 
dominated region. 

This approach also allows a simple theory developed by Stringari |Q for 
linearized collective modes when we use the TF approximation. Taking the 
time-derivative of (1.66|) gives (vco 



0) 



= - V- 



\ dt 



(1.73) 



Using ( 1.7C| ), we have 
dt 



— V [Vex{v) +gncoir) + gSndr^t)] 
m 

--V<5ne(r,t). 

m 



(1.74) 



Combining this last result with ( 1.73| ), we obtain the very useful Stringari 
equation of motion |4G] 



= V 



V5t 



(1.75) 



This describes the collective oscillations of the condensate in terms of a 
single differential equation. As one example, the breathing mode of the 
condensate has a frequency Tito = ^/Zfiujo. This example points out that in 
the TF limit (large Nc), the frequencies are independent of the interaction 
strength and the size of the condensate Nc- 

We also note that using (1.74), (1.75) can be equally well rewritten in 
terms of the superfluid velocity Vc(r, t) defined in ( 1.6^ ). This emphasizes 
that the condensate fiuctuations are directly related to the existence of 
phase fluctuations. Their existence may thus be viewed as "evidence" of 
superfluidity, the latter being always a con seque nce of the phase coherence 
of the macroscopic wavefunction given by ( 1.65 ) [p^ . 

The great thing about the collective oscillations of a condensate in a 
trapped gas is you can "see" them. A beautiful example from MIT is shown 
in Fig. 2 of Ref. |l|. As Ketterle has remarked, these condensates are robust - 
one can kick them, shake them and these "wisps" of Bose-condensed matter 
keep their integrity. 



2.4 Interference of coherent matter waves 

In the pioneering matter wave interference experiments done at MIT using 
a de-tuned cigar-shaped trap |^ , one first destroys the condensate at the 
centre using laser beam. Then the confining trap is turned off and the 
two condensates are allowed to expand and interfere. One observes nice 
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interference fringes at the mid-point, as expected. Using 



^system 

the density is given by 

nir,t) = 1$ 



$A(r,i) + $B(r,i) 



system I 



Na + Nb + 2^/NaNb cos AO, 



(1.76) 



(1.77) 



where A6 = 6a — ds- In the region of interference, the density is low and 
hence interaction effects are small (ie, the gric term is small in the GP 
equation). Asymptotically, the solution of GP equation gives 0{r,t) = ^j-. 
This implies H 



A9 



(- + 1) -(--f)^ 



2ht 



2ht 



mzd 
fit 



2tiz 
WY 



(1.78) 



where X{t) = 2Trht/md. This wavelength is in good agreement with exper- 
imental observations (See Fig. 2 of Ref . ) . 

A condensate described by <&(r,t) may be viewed as a "classical" matter 
wave, as recently emphasized by Pitaevskii and Stringari ji^ . This is quite 
different from ordinary quantum deBroglie waves, since one can ignore 
quantum fluctuations (large Nc). It is also quite different than ordinary 
(classical) macroscopic objects and electromagnetic waves, since $(r,t) is 
described by the GP equation (1.37) which involves Planck's constant h. 
Thus ^{y, t) is a classical object which is described by a quantum equation!! 
This promises to be a challenge for the quantum theory of measurement. 

One can describe a two-component Bose gas (see Section |l|) using coupled 
equations: 



ih. 



5£i 
dt 



2m 
2m 



Vi(r) + ,gi|$i|2 + gi2|$2 



V'2(r) +g2|*2p + .gi2|$i 



$1 

$2. (1.79) 



There are two coupled GP equations for two macroscopic wave functions. 
Extensive studies have been made at JILA using the two atomic hy- 
perfine states of ^''Rb: 



IF = l,mF 



-1 >, \F ^2,mF = l> 



(1.80) 



In particular, one can study interesting interference effects between these 
coupled wavefunctions. Recent work at JILA has used such two-component 
Bose fluids to produce the long sought-for vortex state in one of the com- 
ponents 0. 
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COUPLED DYNAMICS OF THE 
CONDENSATE AND NON-CONDENSATE 



In this Section, we switch our attention from T = (ic, T < 0.4 Tbec) 
to finite temperatures, where iVc and N are comparable in size. We first 
consider how the GP equation of motion for <i>(r,t) is modified. As a first 
step, we could use P, E5|, 



a$(r,t) 
dt 



2m 



+ Ve-^ir) + gn^{r, t) + 2gfi{v, t) 



$(r,i). (1.81) 



The last term takes into account that the condensate moves in the dynamic 
Hartree-Fock (HF) field produced by non-condensate atoms. Immediately, 
one sees this generalized GP equation is no longer closed. It requires a 
theory of the non-condensate fluctuations, as described by n(r, t) = no(r)-|- 
(5n(r, t). 

A simpler version of (|l.81 ) is to treat the effect of the non-condensate as 



a static HF field |38 



3f (pTsi 



2gn{r,t) ~ 2.gno(r). 



(1.82) 



1. This corresponds to treating the condensate moving in a static HF 
field of the non-condensate. 

2. no(r) can be calculated (self-consistently) using the fluctuations of 
<&, as discussed in Section ^. One finds that the depletion of the 
condensate at T = is only a few percent. 

3. This procedure gives reasonable results for the thermodynamic prop- 
erties at finite temperatures, as discussed in the recent literature Q. 
Within the Thomas-Fermi approximation (good for N > 10^ atoms). 



the linearized version of (1.81) using (1.82) leads to the same Stringari 



equation at finite T as the T = result in (1.75). Since the solutions 



of ( 1.75 ) do not depend on the magnitude of the condensate, one con- 
cludes that the collective modes of the condensate will show no tem- 
perature dependence, even though the condensate is being thermally 
depleted. This prediction does not appear to agree with experimental 
results when T > 0.6 Tbec- This suggests that the dynamics of the 
non-condensate has to be included. 



We now go on to determining h(r, t) directly by deriving a quantum 
Boltzmann equation for the single-particle distribution function of excited 
atoms /(p,r,f) and then use: 



n(r,t) 



dp 



/(P,r,t). 



(1.83) 
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This procedure generalizes the approach of Boltzmann (1880's) for a clas- 
sical gas, including the effect of binary collisions. Clearly one must make 
some approximations! One wants, initially, to find a useful kinetic equa- 
tion that builds in just enough physics. Here I will discuss such a quantum 
Boltzmann equation for a trapped Bose-condensed gas at finite tempera- 
tures, which has been extensively discussed by Zaremba, Nikuni and the 
author It is only valid in the so-called semi-classical limit, where it 

is sufficient to work with /(p, r, t). The conditions are 

ksT > gn , ksT > hcuo. (1.84) 

In this domain, one also can assume that the important thermal excitaitons 
can be approximated by simple Hartree-Fock particle-like spectrum: 

2 2 

£(r, t) = ^ + 2g K(r, t) + n(r, t)] + Kx(r) = ^ + U{v, t). (1.85) 

Clearly the resulting kinetic equation is not valid at very low tempera- 
tures, where the thermal excitations are described by a Bogoliubov-type 
spectrum. 

We simply write down our quantum kinetic equation ||^, 

9/(p,r,t) ^ P .v,/(p,r,i)-V,.C/(r,i).Vp/(p,r,i) 

m 



dt 



C22[/]+Ci2[/]. (1.86) 



The right hand side describes how binary collisions effect the value of the 
single-particle distribution function /(p,r,t). The effect of collisions be- 
tween excited atoms in the non-condensate is described by: 

C22[/] = {2-^f%' j / / o'P4'^(P + P2-P3-P4) 

X 5 {£p + Sp^ — — £p^) 

X [(l + /)(l + /2)/3/4-//2(l + /3)(l + /4)]. (1.87) 

This collision integral was discussed in detail in 1933 by Uehling and Uh- 
lenbeck for T > Tbec @1- We recall that creating a Boson gives a factor 
(1 -I- /) and destroying a Boson gives /. In the classical high temperature 
limit, / <C 1 and the collision integral C22 considerably simplifies. 

Where does /(p,r,t) come from in a microscopic derivation of (1.86)? 
Basically we calculate the non-equilibrium real-time single-particle Green's 
functions of the non-condensate field operators (using the Kadanoff-Baym 
formalism |^). This gives (schematically) 

51(1,1') ^(V>(l)Vi(l')), (1.88) 
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where 1 = ri, ii; 1' = r'i,t[. We then express this as 51 (r, t; R, T), where 
the relative and centre of mass coordmates are 

r = ri - r^; R- = ^ (ri +r'i) 

t = ti^t[ r=i(ii+i'i). (1.89) 

Finally we Fourier transform gi (r, t; R, T) to find gi (p, uj;R,T), which gives 
the number of atoms at R, T with momentum p and energy Hlu. The single- 
particle Wigner distribution function is given by 

/oo 
diugi{p,u:-'R,T). (1.90) 
-00 

The Wigner distribution function /(p, R, T) is the quantum generalization 
the classical single-particle distribution function |^ . These remarks should 
indicate how we can go from an equation of motion for the single-particle 
Green's function (within a given self-energy approximation) to a kinetic 
equation for /(p,R, T). We refer to the classic account given (for non- 
Bose-condensed gases) in the book by Kadanoff and Baym for further 
details. This powerful approach was generalized to uniform Bose-condensed 
gases by Kane and Kadanoff [Q, and has been extended to trapped gases 
in recent work ||5l|| . 

In addition to C22 collisions, we also have collisions which involve one 
condensate atom: 

Ci2[/] = j^^^ J dpi J dp2 J dpsSimVc + pi - P2 - P3) 

X 5{ec + Epi ~ e.p2 - Eps) x [(5 (p - Pi) - '5(p - P2) - ^(p - Pa)] 

X K(l + /l)/2/3-«c/l(l + /2)(l + /3)]. (1.91) 

Here the condensate atom has 

energy: = + ^^mvl ; = Vex + gn-c + "^gn 
momentum: Pc = mVc (1-92) 

We note the key difference between C12 and C22 collisions: 

• C22 and C12 conserve energy and momentum in collisions. 

• C12 does not (but C22 does) conserve the number of condensate 
atoms. C12 describes how atoms are "kicked" in and out of conden- 
sate. 



It turns out the generalized GP equation (1.81) is also modified by a 
term related to Ci2[/]. This makes sense, since the C12 collisions modify 
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the condensate wavefunction $(r, i). One finds the new GP equation is 
given by (see also Rcf.p2[): 



ih 



5$(r, t) 



dt 



2m 



14:. (r) + gricir, t) + 2gn{r, t) - ii?(r, t) 



where 



dp Ci2[/(p,r,t)] 



(27r)3 2n,{v,t) 



$(r,t), 
(1.93) 

(1.94) 



More precisely, the dissipative iR term in ( 1.93| ) arises from a three field 
correlation function in the exact equation of motion [see (1.33) and ( 1.35| )], 
and is given by Q 



rir Im 



(1.95) 



We have to solve for /(p, r, t) and <i>(r, t), treating Ci2[/] very carefully. 
We see that there will be an exchange of atoms between the n(r, t) and 
nc(r,t) components through the C12 collisions. We can use these coupled 
equations for a variety of problems. In these lectures, we will consider the 
collective oscillations of the combined system composed of condensate and 
non-condensate. It is useful to introduce two regimes to describe collective 
modes in interacting systems H M : 



27r 



I. Collisionless (produced by mean fields) 
ujtr > 1 or T <C Tfl 

II. Hydrodynamic (produced by collisions) 

uoTji <C 1 or T tr, 



where r^j is some appropriate relaxation time. What should we use for trI 
For a classical gas, this is the collision time [|| 



1 

Tr. 



nav 



(1.96) 



where 



a — Sna^ (for Bosc particles); a ~ s-wave scattering length. 



V ~ average velocity of atoms 
n = density of excited atoms. 



knT 
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Even for a Bose-condensed gas, taking tji ^ Tc is a reasonable first estimate 
[ p7| . To get into the interesting hydrodynamic region {iOTn <SC 1), we need 
small tr, ie, a large density ri or a large collision cross-section a (perhaps 
using a Feshbach resonance, as discussed in Section]^). 

Let us look at the kinetic equation ( [1.86[ ), writing it in the schematic 
form: 

Cf^C22[f]+C\2[f]. (1.97) 

In the coUisionless region, we need only solve £/ = 0. In contrast, in 
the hydrodynamic region, the collisions are so strong they produce local 
equilibrium [||. That is, they force / to satisfy C22[/] = 0. The unique 
solution / of this equation is well-known to be given by 

/(P, r, t) = ■ — , (1.98) 

where v„ is the average local velocity and p, is the local chemical potential 
of the thermal atoms. This local equilibrium Bose distribution involves the 
local variables v„,/i and U, all of which depend on (r,t). 



have [see (1.87) 



Why must / have the form in (1.98|)? To satisfy C22[/i] = 0, we must 



(1 + /i)(l + /2)/3/4 - /l/2(l + /3)(1 + /4) = 0, (1.99) 

and this requires that / be given by the Bose distribution. We have used 
the fact that 

f{x)^-^ = -[fi-x) + l] (1.100) 



and that 

Pi + P2 = Pa + P4 



1 



energy and momentum conservation. 



(1.101) 

2 

where £p = 2m + ^('"' aside, using a kinetic equation is the most 

physical way of deriving the equilibrium Bose distribution. The standard 
approach in statistical mechanics texts based on calculating a partition 
function does not bring out the reason why 

Ibf^ a, \ ■ (1.102) 

However, while the fact that / is given by the local equilibrium Bose 



distribution in (L95) ensures that C22[/] = 0, one finds that Ci2[/] 7^ 0. 



More precisely, we find from (1.91) 



(l + /l)/2/3-/l(l + /2)(l + /3) 



-/3[A-Mc-irn(v„-Ve)"] 



1 



(1 + /l)/2/3. 



(1.103) 
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The expression in the square bracket only vanishes if the condensate and 
non-condensate are in diffusive equilibrium, which requires that 



1 / ^2 



(1.104) 



When we perturb the system, this may not be true, ie, the two components 
may be out of diffusive equilibrium. 

We can now derive hydrodynamic equations for non-condensate by taking 
moments of Boltzmann equation, the standard procedure used in classical 
gases . The first moment gives a continuity equation with a source term: 



dp[Cf = Ci2[f]} 



^ = -V.(nv„) + ri2[/] 



(1.105) 



where 



dp 



/(P,r,i) 



(27r)3- 

;3— /(P'l"'*) 



ri2[/] = 

More explicitly, we find 

ri2[/] : 



(27r)3 TO 
dp 



Cl2[f]. 



(1.106) 



-/3[/i-/Jc-5m(v„-Vc)^l 



1] 



j dpi j dp2 J dp3S{mVc + Pi - P2 - Ps) 



X (5(£c + ei -£2-e3)(l + /i)./2/3 



-/3[/i-pc-|m(v„-Ve)^-l] 



T12 ' 



(1.107) 



We note that T12 is a collision time which descri bes the C12 collisions 
between the C and N.C. at oms. Combining ( 1.105 ) with the continuity 
equation which results from (1.93), 



~dt 



-V-(n,v,)-ri2[/], 

we see that the source term cancels out to give 
d{nc + h) 



dt 



(1.108) 



(1.109) 



Thus our theory gives the exact continuity equation for the total local 
density n = nc + fi. 
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Similarly, one finds 

J dpp[£f = C7i2[/]} 



1. 



dt ' 2 

= -VP(r, t) ~ fiVU{v,t) - to(v„ - v,)ri2[/], (1.110) 
where the kinetic pressure is given by 



^, N Til 



The second moment gives 



dp 

(2^ 



(p- mv„)V(P,r,i). 



(1.111) 



1 I 



u 



ri2[/]. (1.112) 



The detailed derivation of these results is n ot imp or tant h ere The 
main thing is that the hydrodynamic equations ( 1.105 ), ( 1.110 ) and ( 1.112 ) 
can be shown to describe the non-condensate in terms of three new "coarse- 
grained" variables: 

n(r, f), v„(r, i) and P{r,t). 

These are coupled to the two additional variables which describe the con- 
densate: 

nc(r, i), Vc(r,t). 



We no te that the two condensate equations of motion given by ( [l.7C| ) and 
(1.108) are always "hydrodynamic" in form. In contrast, it is only in the 
collision-dominated region that the non-condensate dynamics can be de- 
scribed in terms of a few collective variables. We thus have 5 variables and 
5 equations, which form a closed system. Both components exhibit cou- 
pled, coherent collective motions. This is the essence of two-fluid superfluid 
behaviour |l^, a new unexplored frontier in trapped Bose gases. 

What is new about the two-fluid hydrodynamic equations derived above 
is the role of the source term [/] . In a linearized theory expanded around 
the static equilibrium Bose distribution /o (where ri2[/o] vanishes), one 
finds |, |47) 

ri2[/] - ST,2[f] = -^6,,a^ff, (1.113) 



'12 



where pd. 



Iff 



jl — pLc- Wc find an equation of motion of the kind 

dt T„, 



(1.114) 



where (see Eq.(87) in Ref. ||). 
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Here a involves various static equilibrium thermodynamic functions. The 
new relaxation time (which we can calculate!) determines how fast jl — > 
/Xc, ie, how fast we reach diffusive equihbrium between the condensate and 
non-condensate. We can have 



UJT22 <C 1 
UJT12 <C 1 



required for hydrodynamics (1.116) 



but simultaneously 

LOT^, > 1 (1.117) 

near Tbec, where Uco — > 0. Our hydrodynamic equations predict the exis- 
tence of a new relaxational mode |||, ^ 

lo~-i/t^. (1.118) 

This mode is not included in the standard Landau 2-fluid equations (where 
Ps and pn are assumed to be always in local equilibrium with each other). 

In a uniform gas, the two- fluid hydrodynamic equations give two normal 
mode solutions ^: 

• First sound (oscillation of the non-condensate mainly) 

2 5 Pq kT n 11o^ 

LO = uik , UT ~ - — — ~ . (1.119) 

3 rariQ m 



Second sound (oscillation of the condensate mainly) 



co^U2k,ul^^. (1.120) 
m 

We note the second sound mode is the hydrodynamic version of famous 
T = Bogoliubov phonon mode discussed in Section ^. It is the "soft 
mode" at Tbec- This second sound mode couples to the new relaxational 
mode given in ( 1.118| ) and is damped as a result, the maximum damping 



occuring when ujt^^ — 1. 

In a trapped gas, we can work out the spectrum of hydrodynamic oscilla- 
tions (~ e"*'^*). Both the condensate and non-condensate components have 
the same frequency. The most interesting one is the dipole mode, described 

by 

n(r,i) = no(r-J7„W), JVnW = v„ 

ne(r,i) = n,o(r-J7cW), ^VcW = v,. (1.121) 



One finds there are two modes of this kind |4|, ||] : 
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In-phasc (or Kohn) mode, where ?7„ — rj^ and uj — loq (trap fre- 
quency). It is the finite temperature version of the sloshing mode 



described by (1.63) and (1.64). We note that this mode is generic 
(occuring in both the hydrodynamic and collisionless limit) and is 
not damped Q. 

• Out-of-phase dipole mode, with ?7„ ^ rj^ and in opposite directions. 
The frequency of this mode is different from the trap frequency. This 
out-of-phase mode is of special interest since it is the analogue of the 
out-of-phase second sound mode in supcrfluid "^He. 

We conclude this Section with some remarks: 

1. The specific calculation sketched above is also of interest in the gen- 
eral field of non-equilibrium statistical physics. It describes the de- 
tailed dynamics of a system with a two-component order parameter 
self-consistently coupled to a gas of excitations based on a fully mi- 
croscopic theory. 

2. More work is needed to extend our analysis to low but finite tem- 
peratures and also into the critical region very close to Tbec- In 
both cases, our simple Hartree-Fock particle-like thermal excitation 
spectrum ( |1.85| ) is no longer valid. 

3. The classical kinetic theory of gases has been a rich subject in math- 
ematical physics in the twentieth century, with well-known contri- 
butions by people like Boltzmann, Hilbert, Enskog, Chapman, Uh- 
lenbeck and Burnett. These new equations of motion for a Bose- 
condensed gas promise to yield a lot of new physics in the next cen- 
tury - and surprises, as our work in this Section has already shown. 
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